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ABSTRACT 

It has been recognized long ago that the presence of hot plasma in the inner 
accretion disks around black holes could lead to the neutron production via 
dissociation of helium nuclei. We show that, for a broad range of accretion 
parameters, neutrons effectively decouple from protons and pile up in the 
inner disk leading to the formation of self-sustained halo. This means that 
new neutrons in the halo are supplied mainly by the splitting of helium nuclei 
in their collisions with existing neutrons. Once formed, such a halo can exist 
even if the proton temperature is much lower than the energy threshold of 
helium dissociation. We show that neutron haloes can be the natural source 
of relativistic electrons and positrons, providing characteristic comptonization 
spectra and hard spectral tails observed in many black hole candidates, and 
also giving rise to relativistic outflows. Deuterium gamma-ray line at 2.2 MeV 
resulting from neutron capture is also expected at a level detectable by future 
INTEGRAL mission. Furthermore, the presence of a neutron halo strongly 
affects the dynamics of accretion and leads to the rich variety of transient 
dynamical regimes. 

Subject headings: Accretion, accretion disks - Nuclear reactions, 
nucleosynthesis, abundances - X-rays: binaries 
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1. Introduction 

There are many theoretical and observational indications that the hot plasma of ion 
temperature Tj exceeding several MeV can exist in the inner accretion disks around black 
holes and, possibly, neutron stars. Formation of a two-temperature region with Tj ^> T e in 



accretion disks has been predicted as long ago as in 1976 (Shapiro et al. |1976| ). When the 
energy transfer from ions to electrons is slow as compared with the accretion timescale, an 
accretion flow becomes nearly adiabatic, and the ion temperature can reach virial values. 
This idea has led to the prediction of advection-dominated accretion flows which seems to 



find support in observations of several X-ray transients (Narayan et al. |1996|, |1997|) and 



other sources showing low radiation efficiency and high temperature of the flow. 

Strong observational support for the existence of high-temperature plasma comes 
from gamma-ray and radio astronomy. Very hard spectra extending to MeV region were 
detected from several Galactic X-ray sources presumably containing black holes. Transient 
line-like features around 0.5 MeV indicate an efficient mechanism of positron production 



(Bouchet et al. |1991| , Churazov et al. |1993| ). Non-thermal radio emission and episodic 
ejection of relativistic jets observed in microquasars (Mirabel and Rodriguez |1998| ) are 
definite signatures of ultrarelativistic electrons and/or electron-positron (e~e + ) plasma and 
were interpreted several times as radiation-pressure dominated e~e + outflows (Li and Liang 
TOj BelyaninPOTQ. 



The presence of hot ions in accretion disks leads to the rich variety of inelastic nuclear 
interactions. The most prominent of them is neutron production which primarily occurs 
through dissociation of helium nuclei. This process has potentially important observational 
implications arising e.g. from possibility of radiative neutron capture with subsequent 
formation of gamma-ray lines (Aharonyan and Sunyaev |1984| ). 



The main source of free neutrons in accretion disks is dissociation of helium nuclei 
due to inelastic collisions with energetic ions contained in the disk or free neutrons passing 
through it. Neutron production in accretion disks was considered several times (Aharonyan 
and Sunyaev |1984| Guessoum and Kazanas |1989|) under the assumption that neutrons are 



quickly thermalized in the disk, advected by the bulk plasma motion and transported into 
a compact object on the disk viscous timescale. In this case the neutron fraction in the 
accretion flow is inevitably small and the radiative-capture gamma-ray lines are weak and 
smeared out by a hot rapidly rotating disk. 

However, if plasma in the accretion disk is sufficiently rarefied, then the neutron and 
ion components in the disk effectively decouple, allowing for the accumulation of neutrons. 
This can lead to the formation of neutron haloes with much larger neutron number than 
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in previously considered cases. It is important that, once formed, such a halo becomes 
self- sustained. This means that, even if the ion temperature in the disk falls below the 
threshold for He dissociation by protons or the proton density becomes too small for 
proton-induced dissociation to be efficient, the neutron production is supported by collisions 
of energetic neutrons from the halo with helium nuclei. 

In this paper we propose and explore the scenario of neutron halo formation and show 
that it can lead to much more pronounced observational and dynamical effects as compared 
with previous predictions. 



2. Neutron pile-up and formation of halo 

Let us consider bulk radial motion of neutrons which is the result of angular momentum 
losses caused by elastic collisions with ions (mostly protons) in the accretion disk. One 
directly obtains the following expression: 

V (r) _ d (R g c 2 \ _ c 2 dV n _ R g c 2 

where V n and are the orbital and radial velocities of neutrons, R g the Schwarzschild 
radius of the compact object, and c the speed of light. The time derivative of V n in the 
above expression was substituted by —u pn (y n — Vd)/2, where u pn is the proton-neutron 
collision rate and V d the disk orbital velocity. 

As the orbital velocity of neutrons is larger than that of the accretion disk, the 
excessive centrifugal force has to be balanced by the friction force, which originates from 
the difference of radial velocities of the neutrons and the disk. This gives another relation: 



V 2 - V} 



f (K (r) - vP) . (2) 



R 2 

Here vj r ^ is the radial velocity of the disk and R the current orbital radius. Assuming 
V n ~ Vd — V (Vq is the Keplerian velocity), one finds the expression for the radial velocity 
of neutrons: 

The value (R/V )i , pn in the above equation is the number of collisions a neutron undergoes 
as it completes one orbital revolution, divided by 2n. 

The collision rate v, pn may be expressed in terms of the disk optical depth for 
proton-neutron collisions, r<2, which is equal to the average number of collisions experienced 
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by a neutron passing through accretion disk with Keplerian velocity directed perpendicular 
to the disk plane. Given the heights of ion and neutron "disks", h d and h n respectively, and 
assuming that the cross-section of proton-neutron collisions is inversely proportional to their 
relative velocity, one has v pn = (crV) pn n p (h d /h n ) = (V /h n )((aV) pn /V )n p h d = (V Q /2h n )T d . 
Finally, the expression for the radial velocity of neutrons takes the following form: 

V (r) = {R/2h n ) 2 rj (r) 

4 + {R/2h n frj d ■ 1 1 

Neutrons start to pile up in the inner parts of accretion disk when is considerably 
smaller than vj r \ i.e., when r d ^ 4h n /R. In the limiting case of neutron halo (2h n ~ R) 
the last condition gives r d ^ 1. 



3. Criteria of neutron halo existence 

The neutrons in accretion disk pile up only when the accretion rate M is not too large, 
namely 

Here we used the usual assumption V d ^ = a(h d / R) 2 Vo, and the continuity equation 
to find the proton number density. The ratio h d /R is treated as a parameter. The 
limiting accretion rate M max approximately equals to that corresponding to the Eddington 
luminosity, provided a — 0.1 and h d ~ R. 

With decreasing accretion rate the neutron-to-proton ratio in the disk grows as M~ 2 
and soon exceeds unity, so that collisions with neutrons become the main cause of helium 
dissociation. Actually, it happens at a lower neutron-to-proton ratio, because the neutron 
component tends to have higher temperature than the ion one. Even in the absence of 
additional mechanisms that act specifically to heat neutrons (e.g., precession of neutron's 
orbit in the gravimagnetic field of a rotating compact object), their temperature exceeds 
the ion temperature by a factor 7/5. 

In view of the above-mentioned, the most interesting is the self- sustained mode of 
neutron halo formation, when new neutrons in the halo are supplied solely by the splitting 
of helium nuclei in their collisions with existing neutrons. This mode implies that at least 
one of three neutrons breaks helium nucleus before decaying or being advected into the 
black hole. These two conditions provide, respectively, lower and upper bounds on the 
range of accretion rates in which the self-sustained neutron halo may exist. 
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Let us consider the upper boundary first. To do this, we calculate the average number 
of collisions with protons experienced by each neutron on its way into the black hole: 



//„„/.' ■2;\ll„( :1 ( . \R H 



= -_t^z = _^jr^_ 4 + 

K (r) (*V) pn M\ 



t t V Vn 



-J^ m P , ( 6 ) 



assuming zero helium mass fraction r]H e . In the following considerations we neglect the term 
{R/Vo)v pn as it is important only if M ~ M max . The self-sustained neutron halo requires 



3r] He a d 2a a pn \h d J 

where a pn the proton-neutron scattering cross-section, having the value a pn ~ 700 mb at 
the energies typical for the edge of neutron halo, and a d ~ 100 mb is the dissociation 
cross-section in n — a collisions which, in effect, is the sum of cross-sections for all possible 
spallation processes. It is clear that the self-sustained mode maintains up to M ~ 0AM max 
or even higher; at still higher accretion rates inelastic p — a collisions are frequent enough 
to ensure complete dissociation of helium. 

At very low values of M the neutrons may be advected so slowly that their finite 
lifetime, t n ~ 900 s, comes into play. The decay losses are inevitable (and irreversible), 
but they become a limiting factor only if there are no other processes to compete with, 
for example, evaporation of neutrons from the disk. This means that the ion temperature 
should be considerably lower than the virial value or, in other words, h d R- Keeping in 
mind this restriction, one may derive the absolute lower limit on the accretion rate (the one 
that still permits the self-sustained neutron halo) from the condition 

1 &d 1 

~;VHe Vpntn > o' ( 8 ) 



which gives 



• • 167TCV, Rh n v} r) 4 a pn R f R\ 3/ \ , 

M > M min = v — -^-2 m p = — -—^ 9 - — M max . (9) 

3 a d r]He {o~V) vn t n 3V2 VHeCd ct n \R g J 



The region occupied by the neutron halo is limited in extent by the requirement 
that virial temperature at the edge of the halo approaches helium dissociation threshold 
~ 20 MeV. The latter means R/R g ^ 30. Substituting this ratio in Eq. (§) one finds 
M min ~ 10~ 3 M max for a black hole of ten solar masses and r\n e = 0.2. Thus, there is a large 
parameter space in which the neutron halo may exist. It should be noted that the halo 
can be formed even if M < M m i n , but with a smaller radius. In principle, this marginal 
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solution could exist up to the accretion rates as low as ~ 0.03M m j n , but in this case it faces 
a serious threat from rapidly increasing neutron losses in the immediate vicinity of the last 
stable orbit. This regime, as well as the dynamics of neutron halo in the limit of prevailing 
neutron-neutron collisions, will be considered elsewhere. 



4. Virialization of neutrons in the halo 

As was mentioned above, precession of neutron's orbit caused by the gravimagnetic 
field of a rotating black hole results in additional heating of neutrons and may lead to 
complete virialization of the neutron halo if the precession frequency is higher than or of 
the order of one half of the neutron-proton collision rate, i.e., 

ac '£\\Zn t (10) 



2nR g V R J 2 

where a is the black hole angular momentum in units MR g c/2. Substituting v pn in the 
above equation it is easy to find the corresponding inequality for the accretion rate: 

M max . (11) 



\/27r V R 

Even for a rapidly rotating (a — ► 1) black hole Eq. (|TTD requires mass accretion rate hardly 
consistent with the lower limit M m i n . Therefore, we conclude that the mechanism under 
consideration is unlikely to cause virialization of the neutron halo close to its outer edge. 
However, the virialization is possible in the inner part of the halo at the distance 

( a Mmax \ R . (12) 



5. Long-term variability 

As the neutron halo close to the black hole is likely to be virialized, it produces a 
neutron wind which is nearly isotropic, expands with a velocity comparable to the velocity 
of light, and may carry about a half of all neutrons extracted from helium in an accretion 
disk. Note that helium nuclei completely dissociate when the accretion rate is capable to 
support a self-sustained neutron halo or is comparable with M max , i.e. at any M ^ M min . 
The neutron wind exerts a dynamical pressure on the infalling material via the protons 
originating from the neutron decay. At a distance smaller than t n V w (V w is a typical velocity 
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of neutrons in the wind) this pressure may be expressed in the following form 

p - * -f^W,; < 13 > 

In the case of quasi-spherical accretion (e.g., from stellar wind emanating from a massive 
companion star) the dynamical pressure of the neutron wind prevents the surrounding 
matter from infall provided the dynamical pressure of an accreted plasma at the sonic point 
R s ~ (c/Voo) 2 Rg (Voo is the velocity of the black hole with respect to the stellar wind is less 
than P w . This gives the condition 



As soon as the condition (|I4]) is satisfied, the accretion stops, but the source remains 
active thanks to the matter accumulated in an accretion disk. The disk becomes exhausted 
in a time approximately equal to its viscous timescale at the outer radius, and after that 
the source switches off allowing for a new cycle of accretion from a stellar wind. The quiet 
period lasts for roughly the same time as the active one, until the matter initially deposited 
at the edge of the disk reaches the neutron halo boundary. Such a mechanism results in a 
long-term variability with a period very sensitive to the disk outer size, i.e. to the angular 
momentum per unit mass of the infalling material. The maximum period may be hundreds 
or thousands years. 



6. Proton-neutron bremsstrahlung 

Some of the proton-neutron collisions result in the emission of gamma-quanta, either 
in deuterium formation process or due to bremsstrahlung, first suggested to be an efficient 
source of radiation in proton- neutron plasma by Aharonyan and Sunyaev ( |1984j ). Let us 
define <r 7 so that the energy radiated in each collision of a neutron with a proton constitutes 
on average 0" 7 /<T pn fraction of their relative kinetic energy. The typical value of this ratio is 
of ~ 10~ 4 in our case. Therefore, the energy radiated by a neutron per unit radius is 

dit Opn (Lit 

Here A ~ 2.2 MeV is the binding energy of deuteron, and T is an average temperature 
of protons and neutrons. The radiative efficiency of a disk due to proton-neutron 
bremsstrahlung, as follows from the above equation, is equal to 

(Rjm^ - + A) w^m- (16) 
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This efficiency approaches unity when 



Ai £ 87rr/ He " -m p = £- — M maa . , (17) 



o P n (oV) pn R a a pn \h d J 



assuming A <C 3T and 3T ~ (R g /2R)m p c 2 . At an accretion rate exceeding the above limit, 
luminosity due to proton-neutron bremsstrahlung remains constant independently on the 
accretion rate. 



7. Pion production and hard gamma-ray luminosity of the halo 

As was mentioned above, the inner part of the neutron halo is most likely to be 
virialized. Therefore, neutron-proton and neutron-neutron collisions within R ~ 3R g lead to 
an efficient production of pions. If we assume for simplicity that the mean pion-production 
rate £ n = (a n V), averaged over the velocity distribution of neutrons, is the same for both 
kinds of collisions, the total luminosity in pions is 



1 h 
L n ~ ^{-n 2 n + —n n n p )RlE 7T , 



where R n is the radius within which the pions are produced, E n ~ 140 MeV is the pion 
energy If He dissociation is complete, the density of neutrons in the inner halo can be 
found from the following balance condition: 



r) He M n. 



2,n p (19 > 

where the rate of neutron losses £i oss is mainly determined by the scattering into a black 
hole. Using the expression for n n from Eq. (|19D, the expression for n p derived from the mass 
continuity equation, and assuming h n = R, we obtain pion luminosities associated with 
neutron-neutron and neutron-proton collisions: 

L nn = ^f E ^ (20) 

Sloss ATflp 

Lpn 2vra6ossUJ 1, R 9 cPml ) ' W 
The upper limit on luminosities can be obtained by putting M = M max . This gives 

L- ax ^ 7 x 10 38 ^ He M 10 ( ^) erg/s, (22) 
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L max „ 1Q 38 Miov ^ ^ erg/g? (23) 

where M 10 is the black-hole mass in units of 10 solar masses. Approximately 1/2 of the 
total pion luminosity goes to gamma-quanta and e~e + pairs and can be transformed to 
radiation via inverse Compton scattering of soft X-rays and electromagnetic cascade. 

In the presence of abundant soft X-rays from disk, relativistic e~e + pairs are quickly 
cooled, and the snapshot of electron distribution would demonstrate that the average 
energy is between 50-500 keV. The resulting photon spectrum above 100 keV has a typical 
comptonization pattern. Pair cascade, if developed, leads to a hard power-law tail with 
photon index around 2. Both features are frequently observed in the hard state of accreting 
black-hole candidates. In addition, very hard radiation should be present with photon 
energies up to 70 MeV, but it is expected to contain much less energy than hard X-rays and 
is difficult to observe. 

Neutron halo can be also an efficient source of deuterium line at ~ 2.2 MeV. It 
originates from radiative capture of neutrons that evaporate from the halo and are 
intercepted by the cold outer disk. The total luminosity in the line is (3rju e M A / (2m p ) , 
where (5 is the fraction of neutrons that are captured in the outer disk. Its value depends 
critically on the velocity distribution of neutrons and the size of the halo. For the optimistic 
value (3 ~ 1/2 we have 1,(2.2 MeV) ~ 10" 4 Mc 2 that can be as high as 10 34 erg/s for M 
close to M max - This is detectable by the future INTEGRAL mission from a distance of up 
to several kpc. 



8. Conclusions 

We suggest and analyze in detail the possibility of self-sustained neutron halo formation 
in the vicinity of a disk- accreting black hole (or a neutron star). The neutrons come from 
collisional dissociation of He in the infalling matter. The term "self-sustained" means here 
that, even if the ion temperature in the disk is below the threshold for He dissociation by 
protons, the neutron production is supported by collisions of energetic neutrons from the 
halo with helium nuclei. 

Our analysis shows that halo formation is possible when the mass accretion rate is 
below some critical value. In this case neutrons accrete slower than ions, which leads to the 
neutron pile-up: the density of neutrons can be much higher than proton density. Neutron 
halo exists in a broad range of disk parameters and can lead to the number of observational 
and dynamical effects including self-regulation of accretion and long-term variability of 
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the sources, enhanced hard X-ray and gamma-ray luminosity (also due to proton-neutron 
bremsstrahlung) , 2.2 MeV deuterium gamma-ray line resulting from neutron capture. 
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